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Abstract. A new method is given for the model-space effective interaction. Introducing a new 
operator in place of the Q-box in the Krenciglowa-Kuo (KK) method, we derive a new equation 
for the effective interaction. This equation can be viewed as an extension of the KK method. 
We show that this equation can be solved both in iterative and non-iterative ways. We observe 
that the iteration procedure brings about fast acceleration of convergence compared to the KK 
approach. We also find that the non-iterative calculation reproduces successfully any set of 
the true eigenvalues of the original Hamiltonian. This non-iterative calculation can be made 
regardless of the magnitudes of the overlaps with the model space and the energy differences 
between the unperturbed energy and the eigenvalues to be solved. 



1. Introduction 

In some fields of many-body physics, it is often useful to introduce an effective interaction acting 
in a truncated model space. Much effort has been made both as regards formal theories and 
their applications [1]. 

Among many approaches, we here direct our attention to the Krenciglowa-Kuo (KK)[2] and 
the Lee-Suzuki (LS)[3,4] methods. These two methods are constructed by the use of the so 
called Q-box as the building block of formulation. The KK approach has very simple structure 
and the solution is obtained in an iterative way. If the iteration converges, eigenvalues are given 
for the states with largest overlaps with the chosen model space. On the other hand, the LS 
method reproduces eigenvalues for the states which lie closest to the chosen unperturbed energy. 
The LS method is complicated in structure and higher derivatives of the Q-box with respect to 
starting energy are necessary if one wants to obtain more accurate solutions. 

Both of the two theories yield only certain of the true eigenvalues of the original Hamiltonian. 
This restriction is not, of course, desirable. In a formal point of view, the Q-box itself contains 
information regarding all of the true eigenvalues. For a given model space of dimension d, there 
would be a method of reproducing any d eigenvalues among all the true eigenvalues. 

In this report, we want to show that a method exists for reproducing all the true eigenvalues 
of the original Hamiltonian if the Q-box is given accurately enough. 

2. Krenciglowa-Kuo method 

We assume that the Hamiltonian H is composed of the unperturbed Hamiltonian Hq and the 
perturbation V, i.e., H = Hq + V. The entire Hilbert space is partitioned into the d-dimensional 



model space (P space) and its complement (Q space) with the projection operators P and Q, 
respectively. We note some properties, P + Q = 1, P 2 = P, Q 2 = Q and PQ = QP = 0. We 
here assume that Hq is decoupled between the P and Q spaces as 

H = PH P + QHqQ. (1) 

We further assume that the P-space states have a degenerate unperturbed energy Eq, i.e., 

PH P = E P (2) 

Thus the P-space eigenvalue equation is written with the P-space effective interaction R and 
the P-space eigenstate \<f>k) as 

{EoP + R)\<l> h ) = E h \<j> h ), (fc = l,2,---,d). (3) 

The effective interaction can be constructed by the use of Q-box which is an operator of energy 
variable E and is defined as 

Q(E) = PVP + PVQ E _ 1 QHQ QVP (4) 
The Krenciglowa-Kuo (KK) approach leads to a solution for R as 

d 

Rkk = X)Q(£k)l0k)(0kl 

k=l 
d 

= J2( E k- E o)\<l>k)$k\, (5) 
fc=l 

where (<pk\ is the biorthogonal state of the P-space eigenstate \4>k) in Eq. © defined through 
(0fc|0fc') = ^fefe'- We here note that, as seen in Eq. ([5]), the determination of the effective 
interaction Rkk is equivalent to solving d eigenvalues {E^}. In general, the K-K solution is 
obtained in an iterative way based on Eqs. ([3]) and ([5]). It has been known that if the iteration 
converges, the KK solution yields d eigenvalues of the states with largest P-space overlaps among 
the eigenstates of the original Hamiltonian H. 

We here want to show that there is another way of obtaining the KK solution. We consider 
an eigenvalue equation for a given energy E 

[E + Q(E)]\4, k ) = G k (E)\4> k ), (fc = l 1 2,--. 1 d), (6) 

where {Gk(E)} and {|V>fc)} are the eigenvalues and the eigenstates, respectively. For a given E, 
there are d eigenstates because Q(E) is a (i-dimensional operator. Thus, we have d functions of 
variable E, which we label in order of energy as G\{E) < G2(E) < • • • < Gd(E). 

It may be clear from Eqs. Q and ([6]) that the eigenvalues {E^} in Eq. (|3|) can be given by 
solving equations 

E = G k (E), (k = l,2,---,d). (7) 

We can prove that any solution of Eq. ([7J agrees with one of the true eigenvalues of H and the 
eigenstate \ipk) does with | in Eq. ©. 

The solutions {E k } to Eq. ([7]) are obtained as crossing points of two graphs of y = G k {E) 
and y = E. It should be noted that this procedure of solving eigenvalues {E^} can be made 
independently of the magnitude of the P-space overlap. Therefore we can, in principle, reproduce 
all the true eigenvalues of H. However, as we see in Eq. @ there appear poles in Q(E) when 
E approaches one of the eigenvalues of QHQ. The poles in Q(E) induce also poles in Gk{E). 
Such a situation causes numerical instability in solving Eq. ([7]) with {G k {E)} around the pole 
positions. 



3. Extension of the Krenciglowa-Kuo method 

We define an operator in terms of the Q-box as 



Z(E) = i— - Q(E) -(E- Eq) ■ Qi(E) , (8) 



l-Qi(E) 

where Q\ (E) is the energy derivative of the Q-box given as 

dQ(E) 



Qi(E) 



dE 

PVQ ( E-QHQf QVR (9) 



Hereafter we shall refer to Z{E) as the Z-box. It should be noted that if E = Eq, the Z-box 
agrees with the first-order recursive solution in the LS method[4]. 
The Z-box has the following properties: 

(i) By using Eqs. (gj), (0) and (JSJ), we have 

d ^ 1 

^2Z(E k )\(f> k ){4> k \ = ^2 a /T7 ^kk - Qi(E k ) ■ Rkk \(t>k)(4>k\ 

h=\ k=l 1 ~ Ql{ E k) 

d 

= ^2RKK\4>k)(4>k\ 
k=l 

= Rkk- (10) 

The above fact means that, replacing Q(E) by Z{E) in Eq. ([5]), a new solution for the 
effective interaction R can be derived as 

d 

Rekk = Y,^ E k)\^k){k\ (11) 

k=l 

which we call the extended Krenciglowa-Kuo (EKK) solution. 

In a similar way as in Eq. ([6]) we solve the d-dimensional P-space eigenvalue equation 
for the Z-box for a given variable E as 

{E + Z(E)}\ Xk ) = F k (E)\ Xk ), (k = 1,2, ••• ,d), (12) 

where {F k (E)} are the eigenvalues and {|xfc)} are the eigenstates. We here label F k {E) in 
order of energy as F\{E) < F2(E) < ■ ■ ■ < F^E). In the same way as in Eq. ([7]), the true 
eigenvalues {E k } are given by solving the equations 

E = F k (E), (k = l,2,---,d). (13) 

(ii) We can prove that if E is one of the true eigenvalues {E k } satisfying Eq. (|5|), the energy 
derivative of Z(E) becomes zero, and we have 

This property of the Z-box plays important roles in two ways: 



Figure 1. Dependence of Gj-(E) on E with x = 0.2. The graphs of y = G\{E) and y = G2(E) 
are shown in solid and broken lines, respectively. The direct line denotes the graph of y = E. 

(a) The two equations © and (|12p are not identical. The operator Z(E) is different from 
Q(E) as a P-space operator. Therefore, Eq. (fT3]) based on the Z-box has possibly excessive 
solutions, or spurious solutions, different from the true eigenvalues {E^}. However, we can 
remove easily these spurious solutions according to the condition in Eq. (fT4|) . 

(b) In the same way as in the KK approach we can derive the effective interaction -Rekk m 
an iterative way. The iterative calculation in this approach can be made only by replacing 
Q(E) in the KK method by Z(E). It must be pointed out that, due to the property of 
Eq. ()14p . the iteration based on Z{E) converges quite rapidly as we see in Table 1. The new 
method can be understood as an application of the Newton-Raphson (NR) method which 
is used widely for accelerating convergence in solving non-linear equations iteratively. 

(iii) In the vicinity of the poles in Q(E), Q\(E) has the dominant contribution in Z(E), which 
may be clear from Eqs. dU, ^ and Q. Resultantly Z(E) becomes 

Z(E) « (E-Eo)P, (15) 

and thus there appear no poles in Z(E). Therefore, the functions {Fi t (E)} are continuous 
and differentiable for any E. This fact of the new approach guarantees stability in 
numerically solving the crossing points of two graphs of y = Fk{E) and y = E. We 
want to emphasize here that the procedure of solving {E^} is regardless of the properties 
of the eigenstates of the original Hamiltonian H, i.e., the P-space overlap or the position in 
energy. 

4. A schematic model analysis 

In order to obtain some assessments of the new approach we consider a model problem for which 
exact solutions can be derived easily. The model Hamiltonian we adopt here is given with the 
coupling strength x, the dimension of the entire space n = 4, the dimension of the P space d = 2 
and the unperturbed energy Eq = 1 as 

PHqP = (I 1),pvp=( I £ ) . PVQ = ( t % ) , (16) 



F k (E) 




Figure 2. Dependence of Fk{E) on E with x = 0.2. The graphs of y = F\{E) and y = F 2 (E) 
are shown in solid and broken lines, respectively. The direct line denotes the graph of y = E. 

Table 1. Convergence of the eigenvalue of the lowest-lying state obtained with the KK and 
EKK methods for the model Hamiltonian with x = 0.05. Correct digits with the KK and EKK 
methods are given in n-th iteration step. The starting energies for {E\,E2) are taken to be 
(0.0,0.0). The notation c indicates convergence to more than fifteen decimal places. The exact 
eigenvalue here is 0.8904504858869942. 



n 


correct digits (KK) 


correct digits (EKK) 


1 


0.9- •• 


0.89- • • 


2 


0.89- • • 


0.890450- • • 


3 


0.8904- • • 


0.89045048588699- • • 


4 


0.890450- • • 


c 


5 


0.8904504- • • 


c 


6 


0.89045048- • • 


c 


7 


0.8904504858- • • 


c 



and 

w = (" 5 5 ; X)-^ = ( 3 T 9-5x)- < 17 > 

This model Hamiltonian H was introduced many years ago by Hoffmann et al.[5] and the 
structure of H was investigated in Ref. [6]. We first depict, in Figs. 1 and 2, the dependence of 
the functions Gk{E) and Fk(E), respectively, on the energy variable E with the strength x = 0.2. 
One may see in Fig. 1 that there are two poles in the graph of Gk{E). On the other hand, 
poles disappear in the graph of Fk(E) associated with the Z-box. One observes four crossing 
points between y = Gk(E) and y = E except for pole positions, and six crossing points between 
y = Fk(E) and y = E for k = 1, 2. The spurious solutions in Fig. 2 for y = Fk(E) can be easily 
removed by the condition that the energy derivative at the crossing points should be zero. 



Next we show in Table 1 the results in the iterative method for the KK and EKK solutions. 
The convergence rate in the EKK method is much higher than that in the KK method. The EKK 
method reaches the convergence to fourteen decimal places after only three iterations. However, 
the iterative calculation is not always recommendatory, because convergence conditions are not 
always clear and we cannot control these conditions. In this respect, it would be recommended to 
solve directly the crossing points of y = Gk{E) or y = Fk{E) and y = E. We could reproduce all 
the true eigenvalues of H by calculating the crossing points in a non-iterative way, for example, 
the bisection method. 

5. Conclusion 

In many cases of calculating the effective interaction, iteration and/or recursion methods, for 
example, the KK and LS methods, have been employed. However, one has seen that these 
methods reproduce only certain of the exact solutions. In general, it is quite difficult to control 
convergence of the iteration and/or the recursion. Moreover, these methods do not work well if 
a solution to be solved lies in the vicinity of the pole position. Such a case takes place actually 
in the calculation using the Q-box as in the KK approach. 

In the present report we have proposed a new approach by introducing a new operator, the 
Z-box, and the functions {Fk(E)}. Characteristics of this approach are as follows: (a) The 
eigenvalues in the model space can be obtained by calculating the crossing points of two graphs 
y = Fk(E) and y = E. The functions {Fk(E)} give us all the information on the eigenvalues 
of the original Hamiltonian H. The crossing points can be calculated by using a non-iterative 
method, such as the bisection method. In this case convergence problem does not occur anymore, 
(b) The solutions can be also calculated by applying an iteration method. In this case, the use 
of Fk{E) accelerates convergence compared to the usual KK approach. We want to notice that 
this approach is understood essentially as an application of the Newton-Raphson method for 
solving a non- linear equation, (c) The function Fj t (E) is continuous and differentiable for any 
of E. Instability or difficulty due to the presence of poles does not occur anymore as long as we 
use the function F^E). 

We may conclude that the present formulation on the procedure for obtaining the effective 
interaction removes some difficulties encountered in the usual approaches, such as the KK and 
LS methods. This approach makes us free from divergence of iteration or the pole problem, and 
therefore it may have possibility of applying widely to many-body problems. 
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